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Abstract 

Several transition metals were examined to evaluate their potential for improving the ductil- 
ity of tungsten. The dislocation core structure and Peierls stress and barrier of 1/2(111) screw 
dislocations in binary tungsten-transition metal alloys (Wi_ x TM x ) were investigated using first 
principles electronic structure calculations. The periodic quadrupole approach was applied to 
model the structure of 1/2(111) dislocation. Alloying with transition metals was modeled using 
the virtual crystal approximation and the applicability of this approach was assessed by calcu- 
lating the equilibrium lattice parameter and elastic constants of the tungsten alloys. Reasonable 
agreement was obtained with experimental data and with results obtained from the conventional 
supercell approach. Increasing the concentration of a transition metal from the VIIIA group, i.e. 
the elements in columns headed by Fe, Co and Ni, leads to reduction of the C elastic constant and 
increase of elastic anisotropy k=C^jC . Alloying W with a group VIIIA transition metal changes 
the structure of the dislocation core from symmetric to asymmetric, similar to results obtained for 
Wi-^Re^ alloys in the earlier work of Romaner et al (Phys. Rev. Lett. 104, 195503 (2010)). In 
addition to a change in the core symmetry, the values of the Peierls stress and barrier are reduced. 
The latter effect could lead to increased ductility in a tungsten-based alloy. Our results demon- 
strate that alloying with any of the transition metals from the VIIIA group should have similar 
effect as alloying with Re. 
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I. INTRODUCTION 



Tungsten is the prime candidate for use in the divertor of future fusion reactors because 
of its high melting temperature and resistance to sputtering However, its lack of 

ductility is an impediment to its use. The low-temperature brittleness is a common problem 
for all metals from the VIA group, such as chromium, molybdenum and tungsten 3[. The 
ductility of these metals can be improved by alloying with rhenium, leading to the so-called 
"Re effect" 3^5]. However, Re, is a very rare and expensive element. It is therefore desirable 
to find alternate elements which provide a similar increase in ductility at lower cost. An 
experimental investigation of the range of candidates in the periodic table would be rather 
expensive, but computational materials science methods based on accurate first-principles 
calculations provide a very promising way to narrow the range of possible candidates. 

Several possible mechanisms for the Re-effect have been discussed Sill], and the 
following two mechanisms are selected as most promising in application to monocrystals, i) 
solid solution softening in which an impurity improves mobility of 1/2(111) screw dislocation 
and ii) enhancement of cross-slip in which an impurity modifies the dislocation core structure 
making it easier to cross-slip and increasing the number of possible slip planes. Both of 
these are related to the effect of impurities on the dislocation core and may be amenable 
to investigation by first-principles computational methods even though such calculations 
cannot be used to directly estimate the mechanical properties. Therefore, it should be 
possible to define one or more calculable figures of merit that are related to a materials 
elastic properties and potentially to ductility. Among the potential figures of merit are the 
materials individual elastic constants, Poissons ratio 12|, and Peierls stress; and there is 
strong evidence that these parameters are also directly related to the electronic structure of 
particular impurities. 

The experimentally observed strong correlation between the number of valence electrons 



on the solute atom and the degree o: 



factors in solution softening la, 113]. First principles calculation 18 



softening 131413] points to the importance of electronic 



191 ] indicate that the 



ductilizing effect of Re could orig inate from "direct" improvement of the mobility of 1/2(111) 



screw dislocations 



17, 



20] , either by decreasing the Peierls stress, or enhancing cross- 



slip by changing the slip plane from {110} to {112} which would increase the number 



of available slip planes from 6 to 12 



211 ] . First principles density functional (DFT) 
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calculations demonstrated that alloying with Re in a Wi-^Re^ alloy leads to a transition of 
the 1/2(111) screw dislocation structure from the symmetric core to asymmetric core, and to 



a reduction in Peierls stress 



181 ] . Closely related results were obtained from first-principles 



calculation for Mo alloys with 5d transition metals [20(]. The authors [20[ placed solutes in 
a row along the 1/2(111) dislocation core and calculated the change of stiffness associated 
with moving the row along the (111) direction. According to the results, solutes having 
fewer d electrons (Hf and Ta) increase the stiffness, which authors infer strengthens the Mo 
alloy, whereas those having more d electrons (Re, Os and Ir) decrease stiffness, leading to 
softening. 

The goal of the present investigation is to study possible substitutes for Re in alloys that 



will result in a similar ductilizing effect as Re. The properties o 



alloys, were modeled using the virtual crystal approximation 



W-transition metal (TM) 
22j (VCA). The applicability 



of this approach to modeling elastic properties,structural stability and phonon properties 



.^v-r,,,,, ,„„., ,, d fi^, .v^. n,,„ 



as follows. First, a dislocation is placed in the effective media representing the W-TM alloy. 
The properties of the dislocation in this medium are calculated exactly but for a W-TM alloy 
described by a pseudopotentail which is the weights average of the actual pseudopotentials 
of the impurity and host atoms and the averaged number of electrons. A periodic lattice 
is used with the lattice parameter chosen by obtain zero pressure condition. This approach 
corresponds to a zero order contribution to the electronic system energy expansion with 
respect to the difference between real atomic potentials and the virtual atom. Such an 
approach allows separating the so-called "band structure effects", in this particular case 
filling of d-states, from the effect of local modification of the lattice due to W substitution. 



The latter is not possible in a direct super cell calculation. It was demonstrated [18j that 
dislocation properties are sensitive to modifications of this effective media and alloying 
with Re leads to a sizable reduction of Peierls stress and barrier. Within this approach, a 
comparison of elastic constants calculated using the VCA and super-cell methods provides 
a verification of the VCA accuracy as we demonstrate below. 

We show that alloying with transition metals with a higher number of d electrons (VIIIA 
group) reduces the Peierls stress and barrier for a 1/2(111) screw dislocation. Similar to W- 
Re alloys, the dislocation core symmetry is reduced by alloying with TM from group VIIIA. It 
is demonstrated that the scale of Peierls barrier reductions are similar for all Wi-^TMj. alloys 
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with the same ratio of electrons per atom (e/a). Although, the VCA approach describes 
an influence of band-structure on the properties of screw dislocations, it does not describe 
the discrete nature of alloy structure. However, this approach permits an assessment that 
can significantly reduce the range of possible solute candidates suitable for Re substitution. 
This paper is organized as follows. In section [TT] we briefly review computational details used 
in the calculations. Section l"HI Al describes the electronic structure of tungsten and its alloys 
with transition metals from IVA, VA and VIIA groups. This is an almost complete list of 



possible binary solid solutions [25]. The results of modifications to the elastic constants with 
alloying are discussed in section [Til Bl Section UlI Cl presents our prediction of core structure, 
Peierls barrier and Peierls stress for Wi_ x TM x alloys. Finally, the conclusions are given in 
section [TV] 

II. COMPUTATIONAL APPROACHES 

The electronic structure within the generalized gradient approximation (GGA) of density 
functional theory (DFT) was calculated using the QUANTUM ESPRESSO (QE) pack- 
age 261 ] . The calculation was done using a plane- wave basis set and ultrasoft pseudo- 
potentials optimized in the RRKJ scheme j^. We used the Perdew-Wang 28] exchange- 



correlation functional. The Brillouin zone (BZ) summations were carried out over a 
24 x 24 x 24 BZ grid for the system with one unit cell and 16 x 16 x 16 grid for the supercell 
containing 2x2x2 unit cells, with electronic smearing with a width of 0.02 Ry applied 
according to the Methfessel-Paxton method. The plane wave energy cut off of 42 Ry allows 
reaching an accuracy of 0.2 mRy/atom. As a realization of VCA for the pseudo potential 



method, we used the scheme proposed by Ramer and Rappe 



29] . The elastic constants were 



calculated from the total energies obtained for the set of unit cell deformations 30] . 



We use a periodic quadrupolar arrangement for a 1/2(111) screw dislocation 31] in the 



cell with basis vectors bi = 9ui, b 2 = 5u2 and 63 = U3, where ui = [110], u 2 = [112] and 



M3 = 1/2 [111]. An appropriate choice of lattice vectors [31] reduces the quadrupole cell to 
half the size, and we therefore use a cell with basis vectors equal h\ = (ui+U2+U3)/2, h 2 = b 2 
and h 3 = 63 and 135 atoms. This unit cell contains only two dislocations with opposite 
Burgers vectors (Fig. [TJ). It was demonstrated earlier that this cell size is large enough 



to reproduce such characteristics of the dislocations as Peierls stress and barrier 321438] 
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reasonably well if the elastic interaction correction is included. We are interested in how the 
barriers change with solute concentration and therefore the above correction is not included 



since it's the same for all concentrations. 



for the dislocation calculation, the BZ summation 



was carried out over a 1 x 2 x 8 BZ grid 18] and the initial structure was relaxed until the 
forces were smaller than 0.0005 Ry/A. 

III. RESULTS 

A. Electronic structure 

The electronic density of states (DOS) calculated using the VCA are presented in Figure [2] 
(colored online) for pure tungsten (blue solid line) and two tungsten alloys, one with 6.25 % 
Re (red dashed line) and the second one with 6.25 % Zr (green dash-dot line). Zero energy 
on this plot corresponds to occupation of electronic bands by six electrons. Thus, for pure 
tungsten, zero energy corresponds to the Fermi level. The Fermi energy is placed in the 
pseudo gap between bonding and antibonding d-states. As can be seen from Fig. [2], alloying 
with Re or Zr doesn't produce sizable changes in the DOS, at least in VCA approach. The 
substitution of tungsten with 6.25 % Re leads to an upwards shift of the Fermi energy, shown 
by a red vertical line. The area under DOS between zero energy and the W0.9375Re0.0625 Fermi 
energy corresponds to 0.0625 additional electrons, i.e., occupation of each virtual atom by an 
additional 0.0625 electrons. Correspondingly, alloying with atoms that have a lower number 
of electrons, such as Zr, leads to a downward shift of the Fermi energy, shown by vertical 
green line, and each virtual atom is occupied by 0.0625 fewer electrons. 



ectronic structure in tungsten alloys can be described using the 



Such a modification of e 
rigid band approximation |22j. This approach assumes that alloying does not change the 
DOS but shifts the Fermi level, so the occupation of each virtual atom corresponds to the 
average number of electrons per atom (e/a) in the alloy. Thus, within this approximation, 
all alloys with the same (e/a) value have the same properties. The concentration of different 
alloying elements which have the same (e/a) value can be calculated through the simple 
expression 

* = (1) 

Zj TM — Zj W 

, where Z?m is the number of valence electrons of the alloying atom and Zy/ is the number for 
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tungsten. Although this approach has limited accuracy, it can be very useful for qualitatively 
estimating the change in such properties as elastic constants or Peierls barrier. It should 
be mentioned that alloying of tungsten with TM from VILA, and VIIIA groups, within their 
solubility limits, fills bonding d states with additional electrons and increases strength of the 
bonds. It reflects in the increase of cohesive energy and reduction of Wigner-Seitz radius 
for 4d TM with filling of bonding d states and reduction_of cohesive energy and increase of 
Wigner-Seitz radius with filling of antibonding states 



The VCA approach could be quite inaccurate for describing 



the properties of disordered 



substitutional metallic alloys, especially for transition metals [22] and therefore the accuracy 
of this approach should be analyzed in each particular case. In order to verify the accuracy of 
VCA we compared its results with those obtained using a supercell approach. The supercell 
contained 2x2x2 cubic unit cells with one tungsten atom substituted by a solute atom. 
Thus, the solute concentration corresponds to 6.25 %. In Figure [3j the DOS calculated 
using the supercell (blue line) and VCA models (red line) are presented for WisFe. The 
DOS calculated using the VCA approach is very close to the supercell result. Comparisons 
of elastic properties calculated using the two approaches will be discussed below. 

B. Elastic constants 

The modification of elastic properties in W-Re alloys has been widely investigated both 
experimentally and theoretically 



23 



24j . Here we expand these investigations to a wide 
set of tungsten-based alloys. The results of equilibrium lattice constant, a , bulk modulus, 
B, and elastic constants, Cn, C12, C44 and C = (Cn — C±2)/2, calculated by both VCA and 
supercell methods, are presented in Table [B For the elastic constants, the largest difference 
between calculated and experimental values was obtained for C44 and is about 13 % for 
the VCA calculations and slightly larger for 2x2x2 super cell calculations (SC). The 
agreement between VCA and SC results W alloys with Zr, Ta, and Re is very good and 
differences for elastic constants do not exceed 2 %. Hoever this difference increases for 
alloys with solute atoms with a larger number of valence electron. For example, for the 
Fe-family (Fe, Ru and Os) the difference is about 5 %, while for the Co- (Co, Rh and Ir) and 
for Ni-families (Ni, Pd and Pt) the difference increases to 12-15 %. Since we investigate the 
reduction of Peierls stress and barrier for alloys in which the number of valence electrons per 
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atom (e/a) corresponds to 10 % of Re or smaller the elastic properties of the alloys should 
be described with an accuracy similar to Re of 5 % or less. As discussed above, alloying 
of tungsten with transition metals with a higher number of d electrons, metals from the 
VIIA and VIIIA groups, fills bonding d states. According to the general results obtained 



for transition metals 



22 



391 ]. this filling leads to a decrease of the lattice parameter. The 



calculated results for Wi T Re, x alloys reproduce this general tendency (see Table |TT]) and 



agrees with experiment 



401 ] . The same tendency was obtained in tungsten alloyed with 



other TM except Pd and Pt. For these two elements, the lattice parameter calculated in SC 
slightly increases. Another general trend for all tungsten alloys with solute atoms having a 
higher number of d electrons is reduction of the C elastic constant. The calculated results 
in Table [TT] are also in good agreement with experimental values for the W-Re alloy [40] and 
reflect a reduction in the stability of bcc structure. At concentrations higher than x=0.25% 
the bcc Wi_ x Re x transforms to the cr-phase [41|. Similar to previous results [231 ]. the C 
modulus decreases with Re concentration until it becomes negative at 85 % Re. This change 
in the sign of C corresponds to the dynamic loss of stability of the bcc structure. The C 
elastic constant corresponds to the long-wave transversal phonon branch in [££0] directions 
(T [110][££]) and softening of this phonon mode provides a transition path from the bcc to 
dhcp structure (see discussion in [231]). 

According to the rigid band model discussed in the previous section, the parameter which 
determines the elastic properties of Wi-^TM^ alloys is the number of electrons per atom 
e/a = (1 — x)Zw + xZtm (see explanation of EqJT]). In Fig. [4] the elastic constant C and 
elastic anisotropy A = Cu/C are presented as a function of e/a. The C and A values are 
very similar for all TM from group VIIA and VIIIA except for Fe and Ir. However, even for 
these two elements the deviation from the average values is around 7 %. Additional evidence 
for the importance of the e/a parameter follows from the correlation between the width of 
the stable bcc solid solution region of Wi_ x TM x and the number of valence electrons for 
TM elements from groups VIIA and VIIIA 41] . Alloys with solute atoms having a higher 
number of valence electrons have a narrower solid solution region. 



C. Dislocation core structure, Peierls stress and Peierls barrier. 

In bcc structures 1/2(111) screw dislocations have three nonequivalent core configu- 
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rations, that are referred to as "easy", "hard" and "split" [33|, |42|-|44| - stable and two 
metastable, respectively. These configurations can be obtained for a given Burgers vector b 
by placing the core origin in specific positions. Recently, it was demonstrated that in bcc Fe 
the lowest energy is the symmetric "easy" core configuration, then symmetric "hard" core 
and the highest energy is the "split" core configuration 45]. In a contrast to Fe, for the case 
of W we determined that the symmetric "hard" core configuration decays into the "split" 
core. Thus, in the present publication we assume that the dislocation migrates between 
the "easy" and the "split" core configurations. 

Following, Xu and Moriarty 46[ we marked these sites as 1, 4 or 5 for "easy" cores and site 
2 or 3 for "split" cores in Figure El The calculated structure of the "easy" core configuration 
for a 1/2(111) dislocation is shown in Fig. [5] for pure W, Wo.75Reo.25 and Wo.8sFeo.12, where 
the concentrations of transition metals were chosen to give the same (e/a) value. The circles 
in Fig. [5] represent W atoms looking in the (111) direction and the dislocation structure 
is illustrated by differential displacement (DD) maps 42[. On a DD map the displacement 
of atoms in the [111] direction relative to a neighbor is plotted as an arrow connecting the 
atom with its neighbors. The length of an arrow is normalized so that arrows connecting 
two atoms corresponds to displacement by 6/3. Thus, summing arrows in the circuit around 
a dislocation gives the Burgers vector. For tungsten alloys with VIIIA group TM having 
the same (e/a) value, the core configurations looks exactly the same as for the Re and Fe 
alloys. On the contrary, alloying with Ta or Zr produces the same core symmetry as pure 
tungsten. The same result was obtained for W-Ta by Li et al fl9|. The core structure of 
pure W is symmetric, as shown in Figure [5^, i.e., the dislocation expands equally along the 
six (112) directions, similar to results obtained earlier 18]. 

Alloying with Re or any other TM with a higher number of valence electrons leads to a 
change in dislocation core structure from s ym metric to asymmetric. The core is spread out 
in three (112) directions on {110} planes 46]. There are a six possible (112) orientations 
and the core is thus double degenerate. The transition from symmetric to asymmetric core 
can change the dislocation slip plane 18] . The symmetric core dislocations glide uniformly 
on {110} planes; asymmetric ones glide in a zigzag manner {^tJ and the slip plane changes to 
{112}. Using the core site notation in Figure [5j the symmetric core glide path 1 — > 2 — > 4, 
which is shown by the gray band in Figure [5k will be changed to glide path 1 — > 2 — > 4 — > 
3—7-5, shown by the gray band in Figure Eb 48] . In the present work, the transition from 



S 



symmetric to asymmetric core was obtained in the VCA calculations and this result was 



3Q- 



anc 



49|, 



confirmed for the W-Re alloy in a super cell model calculation reported earlier 

Addition of Re or any TM from group VIIIA leads to the change of dislocation slip p 
and decrease in the value of Peierls stress a p fl8 |. Following a widely used technique 
we apply pure shear strain in the (111) direction which results in a stress that has a zx as its 
main component, the influence of a yx can be neglected since a yx /a zx = 0.04 for pure W. The 
corresponding strain is produced by modification of the basis vector hi along u 3 direction 
hi = 1/2-ui + l/2u 2 + (1/2 + e)u-s. Figure [6] shows the dependence of total energy (see inset) 
per dislocation per Burgers vector and shear stress as a function of strain, e. 

For small strains, i. e. in the elastic regime, the energy increases as e 2 . At larger e 
the energy dependence deviates from square dependence and abruptly drops. This drop in 
energy is caused by a jump of the dislocation core to the next stable "easy" core neighboring 
position (1 — > 4). The corresponding stress and strain are considered to be the Peierls 



stress a p and strain e p . Earlier it was demonstrated 



18[ that the cell size of 135 atoms as 



used in the present work is enough to reach convergence in a p values. Alloying with both 
Re and Fe reduces values of e p and a p . For Fe this reduction is even larger than for Re 
although the solute concentrations corresponds to the same e/a value in both cases. The 
absolute value of a p equal to 1.71 GPa for pure tungsten (Table [TTTT) ) is somewhat lower than 



those obtained earlier 



18|, which were 2.09 GPa using QE and 2.37 GPa using VASP. The 



difference is attributed to the use of different types of pseudo-potentials, ultra-soft in our 
case and norm- conserving (QE) and PAW (VASP) in Romaner et al. 18J], and the higher 
accuracy obtained using our larger plane wave energy cut-off of 42 Ry compared to 30 Ry 
(QE calculations) and 16.4 Ry (VASP calculations) [18] . Moreover, the norm- conserving 
pseudopotential would require an even larger cut-off energy than the ultra-soft one. The 
relative reduction of a p is 20 % for Re and 36 % for Fe. A larger reduction of a p in tungsten- 
iron alloys correlates with the larger reduction of C modulus caused by alloying with Fe 
(Figure H] and Table HI])- These results demonstrate that even if the reduction of Peierls 
stress and C values can be qualitatively understood as arising from a filling of <i-states 
in tungsten by additional electrons from group VII and VIII solutes, some other element- 
specific mechanism exists. It should be mentioned that the solubility limit of Fe in W is 
around 2% and in the concentration of Fe in limited to this value the reduction of Peierls 
barrier is about one half of that of the alloy with 5% Fe. For the W-Ru alloy with a solubility 
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limit of 3 %, a p is reduced by 6 %. This result is consistent with the experimentally observed 



reduction of ductile-brittle transition temperature in W-Ru alloy 



50|. 



The second method used in present work to estimate Peierls stress was originally proposed 
by Nabarro pit], and has been discussed in publications by Groger et al. |48|, |52[. Following 
Nabarro, we measure the system energy as a function of the dislocation location as it moves 
from position 1 to position 2 in Figure [5j It is easy to see that both 1 — > 2 — > 4 ({110} slip 
plane) and 1— >• 2 — > 4 — > 3 — > 5 ({112} slip plane) paths contain a set of dislocation core 
jumps between "easy" and "split" core configurations with the same 1 — > 2 barrier. The 
energy dependence obtained is associated with the Peierls barrier and the Peierls stress is 
given by the maximum gradient of this function 

1 f dV 



(J. 



(2) 



p b \dR, 

where b is the Burgers vector, V is energy per unit length of a straight dislocation and R c is 
the distance along the dislocation core path on the move from 1 to 2. R c equals for position 
1 and \/2a/3, where a is the lattice parameter, for the position 2. We use the drag method, 



also called the reaction coordinate method 



46[, to move the dislocation from 1 ("easy") to 



2 ("split") core configuration. According to this method, the ^-coordinates of atoms in two 
columns around both cores in the modeling cell, shown by purple circles in Fig. [5^, are fixed 
so that z = r c z + (1 — r c )zi, where z and z\ correspond to the z-coordinate of the specified 
atoms in the "easy" and "split" core configurations respectively, and r c = R c /(\/2a/3) is 
the so called reaction coordinate. In the initial implementation 46J, the authors fix the 
^-coordinate of atom which has the largest displacement resulting from the movement of the 
dislocation core. However from our experience the atomic relaxation convergence increases 
if we fix the ^-coordinate of any two of the three atoms surrounding site 2 in Fig. Eh- In 
order to avoid any changes in the elastic interaction between dislocations, both cores are 
moved simultaneously. The results obtained for the Peierls barrier changes are presented 
in Figure [7] for pure tungsten and alloys of tungsten with Re, Fe, Ru and Zr. In order 
to demonstrate that the "split" core configuration corresponds to a metastable state, we 
calculated the potential barrier using a small reaction coordinate step for pure tungsten 
around r c = 1. As can be seen from the results presented in the upper insert in Figure [TJ the 
"split" core configuration corresponds to a local minimum in the potential energy. When 
normalized using the values for pure tungsten, the shape of the Peierls barrier for the alloys 
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is shown to be similar to tungsten and quite insensitive to solute type and concentration, 
see lower insert of Figure [7J In all these alloys, the results yield a single barrier. For alloys 
with Os, Ir, Pt, Ru, Rh, Pd, Co and Ni the total energy was calculated for the "easy" core 
(reaction coordinate, r c = 0) and "split" core (r c = 1). Since the shape of Peierls barrier is 
insensitive to type of solute, the difference between these energies provides an estimate of 
the value of the Peierls barrier as shown in Figure [HJ The reduction of the Peierls barrier 
observed in these calculations correlates well with the reduction in Peierls stress estimated in 
the deformation simulation described above. Alloying with Re or any group VIII transition 
metal at a concentration which corresponds to (e/a) = 6.10 reduces the barrier by 25 %. Fe 
reduces the Peierls barrier even more significantly, whereas Zr increases the barrier. This is 
consistent with the effect obtained for Mo alloys with 5d transition metals in which placing 
a row of solute atoms with higher number of valence d electrons along the dislocation core 
reduce stiffness and a row of atoms with lower number increase stiffness [20I . 

The Peierls stress calculated using the direct deformation method, a p , and from the 
Peierls barrier using EqEJ a p are presented in Table HTT1 For all alloys the Peierls stress values 
obtained from the Peierls barrier are lower than those from the direct method. The difference 
is largest, 14 %, for W alloys with 10 % Re, or 2 % Fe. Thus, the a p results, introduction of 2 
% Fe to tungsten reduces the Peierls stress by 12 %. While for the pure W and W alloy with 
5 % of Fe it is below 5 %. Groger and Vitek [52] connected the discre pan cy in the results 
of two methods with the fact that both drag and Nudged Elastic Band 53] (NEB) methods 
garantee that images of the system are distributed uniformly along the minimum energy 
path and do not implay that dislocation position is distributed al ong this path. In order to 
overcome this problem, the authors introduced a modified NEB [52J method which gives a 
Peierls stress that agrees within 8 % of the directly calculated stress. However we believe that 
the 14 % accuracy obtained in our calculation is reasonably good considering that maintains 
the same trend in a p . It should be metioned that problems with direct application of drag 
or NEB methods can be more substentional. According to modeling result obtained in the 
deformed cell the dislocation core structure is slightly modified with increase of deformation 
and then jumps directly to next "easy" core position. The dislocation core is not observed 
in metastable, "split" core position. Thus the agreemnt between data presented in Table [TTT1 
gives a support to application of drag method to calculation of Peierls stress. 

The reduction of Peierls stress/barrier due to alloying with transition metals with a 
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higher number of valence electrons and the increase for lower number of valence electrons 



is supported by existing experimental data. Hardening was observed in W-Zr alloys 54 1 
and softening in W-Re/Ir alloys Very interesting experimental results for Mo alloys 



with group VIIIA 3d transition metals were presented by Hiraoka et al [13|. The authors 
report a decrease of yield strength with addition of any element from group VIIIA, Fe, 
Co, Ni or Ad Pd. The author concluded that this effect could be well understood in the 
terms number of valence electrons of alloying elements, and that the effect of atomic size 
mismatch is secondary and minor. Although the value of yield strength is not directly related 
to ductility, the observed correlation of properties with e/a value supports the validity of 
our calculation using the VCA approach. 



IV. CONCLUSIONS 



The influence of alloying tungsten with transition metal solutes on the elastic properties 
and properties of 1/2(111) screw dislocation based on electronic structure calculations. In 
comparison with a standard supercell method, we demonstrated that the virtual crystal 
approximation gives a fairly good description of such alloying, especially in the case of 4d 
and 5d transition metals. For the case of alloying with transition metals from group VIIIA, 
the modification of elastic constants, Peierls stress and barrier can be understood within the 
rigid band approximation. This means that alloys with TM concentrations leading to the 
same number of electrons per atoms exhibit a similar reduction of C modulus and elastic 
anisotropy A. Together with this modification of elastic moduli the value of the Peierls stress 
and barrier are reduced by alloying with transition metals from group VIIIA. In addition to 
the Peierls barrier reduction, alloying with any metal from this group changes the "easy" 
core dislocation structure from the symmetric to asymmetric configuration. This similarity 
allows us to conclude that the search for an alternative to Re could be reduced to transition 
metals from the group VIIIA. 
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FIG. 1: The [111] projection of dislocation dipole unit cell, where atoms are shown by 
circles and arrows corresponds to differential displacements. 
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FIG. 3: W . 9375^0.0625 DOS per atom calculated using a 2 x 2 x 2 supercell model, blue 
dashed line, and VCA, red solid line. The Fermi energy corresponds to zero. 
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FIG. 4: C and A as a function of number of electrons per atom e/a calculated for 

Wi-JM, using VCA. 
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FIG. 5: The DD map of "easy" core dislocation in Wi^TMj. calculated using the VCA 
approach. Some of arrows are highlighted in green to illustrate a change of core symmetry. 
The sites 1, 4 and 5 correspond to the stable "easy" core dislocation center, while sites 2 
and 3 indicate the metastable "split" core center. The ^-coordinate of atoms in columns 
shown by filled circles are fixed during calculation of Peierls barrier within the reaction 
coordinate method. The symmetric and asymmetric core glide paths are shown by gray 

bands. 
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100 




Strain (e) 

FIG. 6: Shear stress and total energy, AE, per dislocation per Burgers vector b (shown in 
the inset) as a function of strain, e, for pure W (shown by squares), W .9oR-e .io (shown by 
circles), Wo.95Feo.05 (shown by up pointing triangles) and W .98Fe .o2 (shown by down 

pointing triangles). 
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Reaction coordinate 



FIG. 7: Calculated Peierls barrier, where reaction coordinate, r c , equal for "easy" core 
configuration and 1 for "split" core. The curves are normalized in the center insert. The 
detailed structure of the calculated Peierls barrier in pure W around maximum is shown in 

the insert at upper right. 
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TABLE I: Comparison of lattice parameter (a) in A, bulk modulus (B) and elastic 
constants (Cy) in GPa calculated using the VCA and supercell (SC) approaches for 
Wi^TMj.. The concentration x = 0.0625 corresponds to W15TM. 



TM 


Approach 


a 


B 


a 


C44 


Cn 


C12 


pure W 


VCA 


3.1903 


304 


160 


141 


518 


197 




SC 


3.1896 


308 


157 


139 


517 


203 


Zr 


VCA 


3.2095 


287 


153 


126 


491 


185 




SC 


3.2079 


286 


150 


127 


489 


186 


Ta 


VCA 


3.1946 


298 


161 


136 


513 


191 




SC 


3.1963 


308 


157 


135 


517 


203 


Fe 


VCA 


3.1802 


284 


132 


128 


459 


196 




SC 


3.1716 


300 


140 


132 


487 


207 


Co 


VCA 


3.1788 


284 


125 


115 


450 


200 




SC 


3.1725 


297 


129 


131 


469 


211 


Ni 


VCA 


3.1830 


279 


116 


126 


434 


202 




SC 


3.1743 


295 


126 


131 


463 


211 


Ru 


VCA 


3.1772 


308 


142 


142 


498 


213 




SC 


3.1842 


301 


135 


137 


481 


211 


Rh 


VCA 


3.1710 


309 


135 


141 


490 


219 




SC 


3.1862 


298 


124 


133 


462 


215 


Pd 


VCA 


3.1639 


310 


129 


139 


482 


224 




SC 


3.1901 


293 


116 


131 


448 


216 


Re 


VCA 


3.1857 


308 


151 


144 


510 


207 




SC 


3.1859 


308 


150 


143 


508 


208 


Os 


VCA 


3.1809 


312 


143 


146 


503 


216 




SC 


3.1848 


307 


136 


144 


489 


216 


Ir 


VCA 


3.1736 


316 


147 


148 


513 


218 




SC 


3.1865 


304 


125 


140 


470 


221 


Pt 


VCA 


3.1687 


317 


132 


148 


494 


229 




SC 


3.1917 


299 


113 


135 


449 


223 
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TABLE II: Experimental and calculated (based on VCA) values of lattice parameter (a) in 
A, bulk modulus (B) and elastic constants (C^) in GPa for Wi^TM^. The concentration 
for all TM except Re is chosen to keep number of electrons per atom (e/a) = 6.10. 



TM 


X 


a 


B 


a 


C44 


Cu 


C\2 


experiment [4QJ 


W 


0.0000 


3.1659 


314 


164 


163 


533 


205 


calculated, virtual crystal approximation 


Re 


0.0000 


3.1903 


304 


160 


141 


518 


197 




0.0300 


3.1880 


306 


156 


142 


514 


202 




0.0500 


3.1866 


307 


153 


143 


512 


205 




0.1000 


3.1831 


310 


146 


145 


504 


212 




0.3000 


3.1689 


322 


130 


156 


495 


235 




1.0000 


3.1317 


350 


-30 


163 


309 


370 


Zr 


0.0250 


3.1972 


297 


161 


136 


511 


190 


Ta 


0.0500 


3.1936 


299 


162 


137 


515 


191 


Fe 


0.0250 


3.1849 


296 


148 


136 


494 


198 


Ru 


0.0250 


3.1846 


306 


153 


141 


510 


204 


Rh 


0.0167 


3.1846 


306 


153 


141 


509 


204 


Os 


0.0250 


3.1862 


307 


153 


143 


511 


205 


Ir 


0.0167 


3.1849 


308 


154 


143 


512 


205 



TABLE III: In order to compare the reduction of Peierls stress calculated using the direct 
deformation method, a p in GPa, and from the Peierls barrier using Eqf2J a p in W, 
W . 90 Reo.io, Wo.95Feo.05, Wo.98Fe .o2 and Wo.97Ruo.03- 







a p 


w 


1.71 


1.64 


W0.95zr0.05 


2.18 


1.82 


Wo.90Reo.10 


1.37 


1.15 


Wo.95Feo.05 


1.09 


1.05 


Wo.98Feo.02 


1.65 


1.40 


Wo.97Ruo.03 


1.60 


1.35 
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